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ABSTRACT 


Finite element thermal stress analysis was performed on a rectangular titanium honeycomb-core sandwich panel 
which is subjected to thermal load with a temperature gradient across its depth. The distributions of normal stresses 
in the face sheets and the face-sheet/sandwich-core interfacial shear stresses are presented. The thermal buckling of 
the heated face sheet was analyzed by assuming the face sheet to be resting on an elastic foundation representing the 
sandwich core. Thermal buckling curves and thermal buckling load surface are presented for setting the limit for 
temperature gradient across the panel depth. 


INTRODUCTION 


There are a variety of candidate hot-structural panel concepts for application to hypersonic aircraft such as 
the national aero-space plane and the high-speed civil transport. One of the promising candidates is a titanium 
honeycomb-core sandwich panel (both face sheets and the sandwich core are made of titanium material). During 
high-speed flight, the hot side of the panel will be heated to higher temperatures than the cooler side, resulting in 
the buildup of a severe temperature gradient across the panel depth. Under this type of thermal loading, several 
phenomena could happen. The typical ones are (1) debonding between the heated face sheet and the sandwich core 
if the interfacial shear stresses exceed the shear bonding strength, and (2) thermal buckling of the heated face sheet 
if the temperature gradient across the sandwich panel is too severe. 

This report presents thermal stress analysis of the titanium sandwich panel, and examines whether these 
phenomena will take place under the proposed service heating conditions and laboratory test conditions of the 
sandwich panel. 
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length of sandwich panel, in. 

Fourier coefficient of assumed deflection function for w( x, y) 
width of sandwich panel, in. 


face sheet flexural rigidity, 
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modulus of elasticity of face sheet material, lb/in 2 

moduli of elasticity of sandwich core, lb/in 2 

shear modulus of face sheet material, lb/in 2 

shear moduli of sandwich core, lb/in 2 

sandwich-core depth, in. 

number of buckle half waves in x-direction 

panel compressive load in x-direction, lb/in. 

panel compressive buckling load in x-direction, lb/in. 

panel compressive load in y-direction, lb/in. 

number of buckle half waves in y-direction 

structural performance and resizing 

thickness of face sheets, in. 

strain energy of plate bending, in-lb 
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A C/2 

strain eneigy of elastic foundation, in-lb 

AW 

work done by external forces, in-lb 

w 

face sheet out-of-plane deflection, in. 

x,y 

rectangular Cartesian coordinates, in. 

ot c 

coefficient of thermal expansion of sandwich core, in/m.-°F 

Olt 

coefficient of thermal expansion of face sheet, in./ia-°F 

P 

foundation modulus, lb/in.-in 2 

7c 

weight density of sandwich core, lb/in 3 

Is 

weight density of face sheet material, lb/in 3 

Vs 

Poisson ratio of face sheet material 

Vxy i Vy Z , V X z 

Poisson ratios of sandwich core 


normal stress in face sheet in the x direction, lb/in 2 

Tcxz 

interfacial shear stress, lb/in 2 


DESCRIPTION OF PROBLEM 


Figure 1 shows the rectangular honeycomb-core sandwich panel under consideration. The panel is subjected to 
heating with a temperature gradient across its depth. Under such thermal loading, debonding between the heated 
face sheet and the sandwich core or thermal buckling of the heated face sheet could occur if the thermal gradient 
across the sandwich depth is too steep. 

The problem is to perform thermal stress analysis and thermal buckling analysis of the sandwich panel and 
to examine whether these two phenomena might take place under the proposed service heating condition of the 
sandwich panel. 


THERMAL STRESS ANALYSIS 

In the thermal stress analysis, the sandwich material will be considered linearly elastic (neglecting thermal creep 
effect), and the edges of the sandwich panel will be clamped. 

The thermal stress analysis will yield compressive stress distribution in the heated face sheet and will yield the 
interfacial shear stress distribution built up at the interface between the heated face sheet and the sandwich core. The 
face sheet peak compressive stress will then be compared with the thermal buckling stress of the heated face sheet, 
and the peak interfacial shear stress with the interfacial shear bonding strength. 

The structural performance and resizing (SPAR) finite element computer program (ref. 1 ) was used in the thermal 
stress analysis. Figure 2 shows a SPAR finite element model set up for a thin strip of the sandwich panel along the 
y - b/2 line. Each of the face sheets was modeled with only one layer of elements in thickness direction, and the 
sandwich core with six layers of elements in thickness direction. Because of steep stress gradients near the panel 
edges, the element densities were increased in two steps near the panel edge for obtaining smooth stress distributions 
in this region. 


THERMAL BUCKLING ANALYSIS 

When the upper face sheet is heated more than the lower face sheet, the upper face sheet will be under 
compression (upper part of fig. 3). If the sandwich core is considered as a spring, the situation could be approx- 


2 


imated by the buckling of a single face sheet on elastic foundation, as shown in the lower part of figure 3. 

In the thermal buckling analysis of the heated face sheet, the energy method will be used. Let A U \ , A U 2 . A W 
be, respectively, the strain energy of plate bending, strain energy of elastic foundation, and the work done by the 
external forces. At thermal buckling the following energy balance must hold: 

AUi + AUi = AW (1) 


These three energy terms may be expressed in terms of plate deflection w( x,y) as follows: 
Strain energy of bending A U \ : 
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where D is the face-sheet flexural rigidity given by 

D= 12(1 -v 2 ) 

where E Si v Si t s are, respectively, face sheet modulus of elasticity, face sheet Poisson’s ratio, and face sheet thickness. 
Strain energy of elastic foundation (A Ui ): 
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A [/2 = — [ [ w 2 dxdy 
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where /? is the foundation modulus, given by 


p= 


E z 


(4) 
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where E z and h c are, respectively, modulus of elasticity and depth of the sandwich core. 
Work due to external forces (A W): 
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where N x and N v are, respectively, the compressive panel loads in the x- and y -direction in the face sheet. 

For simply supported edges, the face sheet deflection surface may be represented by the following double Fourier 
series (ref. 2): 


u>(x,y) = 52 “mn sin sin 


mry 


(7) 


where a mn is the deflection coefficient, a and b are the length and the width of the sandwich panel, respectively, and 
m and n are the number of buckle half waves in the x- and y -directions, respectively. 

Substituting equation (7) into equations (2), (4), and (6) yields 
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In the light of equations (8) through (10), the buckling equation (eq. 1) takes on the form 
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This expression becomes a minimum if all the coefficients o m », except one, are taken to be zero (ref. 2). 
Thus, the minimum buckling load is given by 
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For thermal loading with constrained edges, N x — N v , and equation (12) becomes 
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For a fixed o and 6, the value of ( N x ) e T will become lowest at particular combinations of m and n. 

To find the functional relationship between m and n which gives lowest value of (N x )a-i we differentiate equa- 
tion (13) with respect to m or n, and set the derivatives to zero: 


from which we obtain 
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which is the equation of ellipse in the m-n space shown in figure 4. 

Substituting equation (15) into equation (13) yields the lowest buckling load of ( N x ). 

(N x )cr = Is/DP 


< CT ■ 


(17) 


Equation (17) is independent of panel aspect ratio, and is subjected to the condition in equation (15) or (16) which 
relate mton 


For a square panel a = b, and equation (15) becomes 


2 2 / P 

m +n = ^D 


( 18 ) 


which describes a circular arc in the m-n space for which m > 1 and n > 1 . 
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NUMERICAL EXAMPLES 


The titanium honeycomb-core sandwich panel, which has been designed and fabricated for testing in the NASA 
Ames Research Center’s Dryden Thermostructures Research Facility, has the following geometry and material prop- 
erties. Geometry: 


o = b = 24 in. 
h c = 0 .69 in. 
t t = 0 .06 in. 


Material properties of titanium face sheets: 


ITEMS 

300°F 

900° F 

EZ lb/in 2 

15.1 x 10* 

11.8 x 10 6 

G„, lb/in 2 

5.7634 x 10 6 

4.5038 x 10 6 


0.31 

0.31 

a s , in/in.-°F 

5.1 x 10- 6 

5.6 x 10~ 6 

lb/in 3 

0.16 

0.16 


Material properties of titanium honeycomb core: 


ITEMS 

600°F 

E x , lb/in z 

2.7778 x 10 4 

lb/in 2 

2.7778 x 10 4 

E x , lb/in 2 

2.7778 x 10 5 

G xy , lb/in 2 

0.00613 

G vz , lb/in 2 

1.81 x 10 5 

G xx , lb/in 2 

0.81967 x 10 5 

Vxy 

0.658 x 10 -2 

V V z 

0.643 x 10" 6 

Vxz 

0.643 x 10~ 6 

a c , in./in.-°F 

5.35 x 10" 6 

7c, lb/in 3 

3.674 x 10" 3 


The sandwich panel will be subjected to the proposed service thermal loading condition of 300°F on the lower 
face (Tt = 300 °F), and 900°F on the upper face (T u = 900 °F). Thus, it will be under the temperature gradient of 
(T u — Ti ) /h c = 600 °F jh c . The four edges of the sandwich panel will be restrained with resisting moments to keep 
the panel flat during thermal loading. This edge condition will induce the highest compressive stresses in the upper 
face sheet (away from the edges). 

In the finite element computations, dissimilar material properties will be used for the lower and the upper face 
sheets. However, for the honeycomb core, the averaged material properties at 600°F will be used. 

RESULTS 

Figure 5 shows the distributions of the magnitude of the normal stress Of x in the face sheets and the interfacial 
shear stress in the y = b/z plane, calculated from the finite element analysis. The normal stress a j x is maximum 
at the panel center with the value (df x ) max = 23,193 lb/in. and remains almost constant up to 6 in. from the panel 
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center, and then monotonicaUy decreases to zero at the panel edge. On the other hand, T a , is zero in the panel center 
region up to 6 in. from the panel center, and increases rapidly to its peak value of (r cxz ) max = 706 lb/in 2 near the 
panel edge. It then decreases sharply to zero at the panel edge. For this low value of interfacial shear stress, it is 
unlikely that the interfacial shear debonding will take place under the proposed service thermal loading. 

Next, we will examine whether the compressive stress generated in the upper face sheet will cause it to buckle. 
The compressive panel load in the upper face sheet is given by 

N X = N V = (c 7f x )rnaxt t = 1392 lb/in. 

From equation (17), the lowest buckling load for the present case is 

( N x )cr = 2 \/~D]3 = 19,452 lb/in. 


for all combinations of m and n. 

Thus, the induced compressive panel load in the upper face sheet is approximately 7 percent of the predicted 
buckling load, and the possibility of face sheet buckling is quite remote. Figure 6 shows ( N x )cr ( e 9- (13)) plotted as 
a function of m for n = 1 and for n = m. For n = 1 the lowest buckling load occurred at m = 49 , and for n = m it 
occurred at m = 35 . The level of compressive panel load induced in the upper face sheet is shown for comparison 
with the predicted buckling loads. Figure 7 shows the buckling load surface plotted in the m-n space. The lowest 
compressive buckling load points form a circular arc at the height of ( N x ) „ = 2 \JD$ measured from the m-n plane 
(horizontal plane). 


CONCLUSIONS 

Thermal stress analysis was performed on a titanium honeycomb-core sandwich panel subjected to heating with a 
temperature gradient across its depth. The magnitude of the interfacial shear stress induced under the proposed ther- 
mal loading was relatively low and would not cause interfacial debonding between the face sheets and the sandwich 
core. 

Also, thermal buckling of the heated face sheet was analyzed by assuming the face sheet to be resting on an 
elastic foundation which is the sandwich core. The calculated compressive thermal stress induced in the heated face 
sheet was approximately 7 percent of the compressive buckling stress predicted from the buckling analysis. Thus, 
the heated face sheet is unlikely to buckle under the proposed thermal loading. 

Ames Research Center 
Dryden Flight Research Facility 
National Aeronautics and Space Administration 
Edwards, California, September 27, 1990 
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Tj = 300° F 
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Figure 2. Finite element model for a thin strip of sandwich panel along y = j plane. 
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Figure 3. Representation of thermal buckling of heated face sheet of sandwich panel 
by compressive buckling of thin plate on an elastic foundation. 



Figure 4. ( N x ) cr = Minimum curve for regtangular sandwich panel. 
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Figure 5. Face sheet normal stress and interfacial shear stress induced 
in titanium sandwich panel under thermal load. 
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Figure 7. Thermal buckling load surface in m-n space for titanium 
honeycomb-core sandwich panel face sheet, a = b. 
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